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Wave force measurements were taken on a horizontal 

30.48cm (12.0") diameter circular cylinder at the Oregon 
State University Wave Research Facility to determine the in- 
fluence of a plane boundary on these forces. The longitudi- 
nal axis of the cylinder was oriented perpendicular to the 
incident gravity wave flow. H/T? ranged from 0.0085 to 


505 ft/sec’. h/t? ranged from 0.078 to 2.57 ft/sec. 


Maximum values of Reynolds numbers ranged from 2.0 x Mo Eo 
Ze x 10°. Force coefficients of inertia, lift and drag for 
a horizontal circular cylinder varied according to the proxi- 
mity of a plane boundary, expressed as e/D, and the relative 
horizontal displacement of the water particles, expressed as 
A/D. Values of Cr increased from Cr = 1.9 to Cr = 3.5 as 

the cylinder approached the plane boundary. Agreement with 
theoretical values for a uniformly accelerated flow was good. 
L varied from ns 0 

ele) Cc, = -5.3 as the cylinder approached the plane boundary. 


For A/D < 1.5 and e/D > 0, values of C 





Agreement with potential flow theory was good. For 1.5 < 


A/D < 5.0 and e/D = 0, Cy. See hOr 5) ar) < 527 and 


e/D > 0, Cy. increased from Cr. = () to Cr = 2.6 as the cylin- 


der approached the plane boundary. For 1.5 < A/D < 5.7 and 
0 < e/D < 1.0, large alternating positive and negative lift 


forces occurred at twice the wave frequency. As the cylin- 


~ 5 


der approached the plane boundary for Re = 1.3 x 10°, Ch in- 


creased from Cy = 0.66 at e/D = 1.0 to Cy = 1.9 at e/D = 
0.17, and then decreased to Ch = 0.93 at e/D = 0. The pos- 
Sibility of a phase difference between the ambient velocity 
at the cylinder and the drag force was shown. Variations 


in water depth produced variations in C, and C,. As h/D 


IL 
varied from h/D = 2.35 to h/D = 7.8, Cr evaluted from verti- 
cal forces increased from Cs 1.8 to Cy = 3.08 and Cy 
evaluated from horizontal forces increased from C, =) 2a Beo 


Cr =EoUOge Or A, D s 1.5, as h/) varied from h/D = 2.3550 


h/D = 7.8, Cy varied from C,. = -~4,29 to C, = f= 526. sah 


A/D > 1.5, water depth had no apparent effect on C, or Cp: 


L 


The advantages of using the Stream function wave theory in 


lieu of the Airy wave theory were demonstrated. 
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WAVE FORCES ON A HORIZONTAL CIRCULAR CYLINDER 
NEAR A PLANE BOUNDARY 


Le LCE ObUeT LON 


In recent years engineering design of ocean structures 
has become more highly developed due to the increased demand 
for ocean resources and offshore storage and living facili- 
ties. Much research is being conducted to better understand 
the dynamics of the ocean, the fluid flow around structures 
placed in the ocean, and the response of the structures to 
the ocean waves and currents. 

To establish a useful foundation for the design of 
ocean structures, at least two areas of research are impor- 
tant. First, research can be conducted in the field to de- 
velop an analytical or numerical model which will predict 
the ocean kinematics and dynamics. The theoretical results 
can be correlated with measured forces on a prototype struc- 
ture. Second, research can be conducted in the laboratory 
to correlate the known or well-controlled flow characteris- 
tics with measured forces on a model. The experimentation 
and conclusions contained herein reflect a laboratory model 
study conducted at the Oregon State University Wave Research 
Facility through the National Oceanic and Atmospheric Admin- 
istration under Sea Grant Contract 04-5-158-2. The forces 
on a horizontal circular cylinder near a plane boundary due 


to gravity waves were specifically considered in this 





Z 
investigation since circular cvlinders are commonly used in 
ocean structures and can aiso represent pipeline networks 
between these structures and shore. 

The study of hydrodynamic forces on cylinders can be 
classified into at least four major flow regimes: (1) steady 
flow at subcritical Reynolds numbers; (2) steady flow at 
Supercritical Reynolds numbers; (3) unsteady flow at sub- 
critical Reynolds numbers; and (4) unsteady flow at super- 
critical Reynolds numbers. Understandably, experimental 
results have been most abundant for the steady flow regions. 
Most laboratory modeling has produced few results at super- 
critical Reynolds numbers. 

Analysts have endeavored to correlate the hydrodynamic 
forces on cylinders in unsteady flow to the kinematics of 
the flow field. The most widely used approach has been the 
Morison equation which represents the total force on a 
structure as the linear sum of an inertia and drag force. 
These component forces are related to the flow kinematics 
through certain dimensionless force coefficients which may 
be determined analytically or experimentally. 

The purpose of this thesis is to report on the results 
of research into the effect of a plane boundary on the hydro- 
dynamic forces on a horizontal circular cylinder whose longi- 
tudinal axis iS oriented perpendicular to incident gravity 


wave flow. The following topics are specifically examined: 
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The contribution of potential flow effects and wake 
effects on the magnitude and direction of the hydro- 
dynamic forces; 

The relation between the force coefficients of in- 
ertia and lift and the proximity of the cylinder to 
a plane boundary; 

The relation between the coefficient of drag and 

the proximity of the cylinder to a plane boundary; 
and 

The effect of variations in water depth on the 


force coefficients of inertia, lift, and drag. 





II. REVIEW 


Past theoretical and experimental investigations of 
both steady and unsteady flow around cylinders will be brief- 
ly reviewed. This review will facilitate the description of 
the mechanics of unsteady flow across a horizontal cylinder 


and of the resulting hydrodynamic forces on the cylinder. 
Wave Force Analysis 


Wave forces on submerged cylinders may be estimated by 
analytical methods or empirical methods. Closed solutions 
of a mathematical model resulting from analytical methods 
are infrequently obtained. These solutions arise from some 
form of the Navier-Stokes equations when viscous effects are 
considered or from the Euler equation when viscous effects 
are ignored. Empirical methods of force analysis are used 
more frequently however since the oscillatory flow around a 
cylinder can be a complex combination of potential flow 
characteristics and the real flow characteristics of free 
stream turbulence and wake formation. The most extensively 
used empirical method has been the so-cailed Morison 
equation. 

Morison, O'Brien, Johnson and Schaaf (8) conducted ini- 
tial research specifically related to wave forces on verti- 
cal piling. Their approach represents the total force on a 


structure as a linear sum of drag and inertia terms. These 





terms are in turn related to the kinematic flow field and 
certain hydrodynamic force coefficients. The Morison eaqua- 
tion shows the total horizontal force per unit length ona 


cylinder to be 


2. 
F=C Se een ULUL (1) 


where 


Dp ST 7 
2 


The coefficients of inertia and drag are in most cases ex- 
perimentally determined. However, the inertia coefficient 
may be analytically determined from potential flow theory. 

The above force coefficients have been computed by cor- 
relating measured forces with theoretical kinematics as in- 
vestigated by Dean and Aagaard (3) using the Stream function 
wave theory or by correlating measured forces with measured 
kinematics as investigated by Keulegan and Carpenter (7) and 
more recently, Sarpkaya (il). 

In order to establish design values of the hydrodynamic 
force coefficients, investigators have sought to find some 
dependence of these coefficients on certain flow parameters. 
For example, for a cylinder in an infinite fluid in steady 


flow, four such parameters are the Reynolds number, cylinder 
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roughness, free stream turbulence, and boundary conditions. 
They alter the flow configuration and hence the drag coeffi- 
Cient. The Reynolds number is a Similarity parameter which 


expresses the ratio of inertia forces to viscous forces as 


Re = OUD | 
u 


A distinct dependence of Ch and the Reynolds number has been 
found for a cylinder in an infinite fluid in steady flow. 


Roshko (9) shows that Figure 1 can be divided into five 


regions: 
1. Subcritical where Cy = lee LOT HO < Re < 2 x tie 
2. Lower transition (2 x 10 < Re < MO ys 
aa, | SUupM@reritical (5 x 10> < Re < ie) where C. = 0.3; 


D 

4, Upper transition (10° < Re < 3.5 x We = and 

See transcratical (Rw, > 3.5 x 10°) where Cy = 0.7. 

In the presence of a plane boundary or when a cylinder 
is immersed in a boundary layer, the coefficient of drag is 
altered. For 0.1 < e/D < 1.7, as defined in Figure 2, 
Wilson and Caldwell (14) found that Ch 
1.8 with Re = 3.3 x 10° and from 1.1 to 1.25 with Re = 


varied from 1.7 to 


preo6 xX 10 For a cylinder resting on the bottom, Beattie, 
Brown and Webb (1) show an approximate value of Cy = 0.4 for 
Reynolds numbers greater than 0.5 x Wo. 

Wilson et al. (14) also show that the coefficient of 


lift varies with the presence of a plane boundary in steady 


flow where the coefficient of lift is defined as 






1 10 10° 102 10% 10° 


Reynolds Number, Re 


Figure 1. Drag coefficient for circular cylinder as a func- 
tion of Reynolds number (from Ref. 9). 
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10 19 


Figure 2. Cylinder cross-section close to plane boundary. 









F 
Cc, = 2. (4) 
4 ODU 
ROE Re = 3.3 x 10°, C. = 0.57 atve/D = 0.02 and C, = 0.33 at 
ype 0.25. For R&= 5.66 x 10°, C. = 0.28 at e/D = 0.02 
and C. = 0.17 at e/D = 0.25. Rouse (10) shows that for Re = 


L 
5 


eo xs 10, C. = 0.50 at e/7D = 0. 


L 
Therefore, for a stationary cylinder placed in steady 
flow, the force coefficients of drag and lift vary with the 
Reynolds number and the proximity of a plane boundary. In 
a uniformly accelerated flow, the coefficient of inertia 
for a cylinder is also a function of the proximity of a 
plane boundary. It is a function of the Reynolds number 
only when a wake forms. The wake alters the effective shape 
of the cylinder which in turn alters the coefficient of 
inertia. 
Within the unsteady flow associated with the passage 
of a wave, the cylinder experiences effects from potential 
flow conditions and real flow conditions. For potential 
flow, the flow field is accelerating and irrotational ex- 
cept in the thin boundary layer at the surface of the cylin- 
der. For real flow, the development of a uniform wake or 
an alternating wake due to vortex shedding alters the flow 
field from potential flow. Force coefficients determined 
under unsteady flow conditions therefore, respond to either 
potential flow effects or wake effects. And as with steady 


flow, it is equally important to determine whether these 






9 
force coefficients vary with the Reynolds number, the proxi- 
mity of a plane boundary, or other flow parameters such as 
is 


UT 
“a= discussed below. Thus, conditions may occur where C 


dependent on wake effects and Cr and C 


D 
are dependent on 


L 
both wake and potential flow effects. Since the presence of 
a plane boundary alters both the characteristics of a wake 
and the potential flow field, then conditions may also occur 
where these force coefficients are dependent on the location 
of the plane boundary. 

Using measured forces on vertical piling from a field 
project correlated with the theoretical kinematics provided 
by the Stream function wave theory, Dean et al. (3) showed 
that within the range where their draq coefficient analysis 


was well-conditioned (Re > 6.0 x Mo? - C. decreased with 


D 
increasing Reynolds numbers. For Re < 3 x We, they found 


that the steady state value of C, = 1.2 from Roshko (9) was 


D 
consistent with data analysis results. Values of Cr were 
shown to be constant at Cy x= 1.5 for all Reynolds numbers 
where the inertia coefficient analysis was well-conditioned. 
In a laboratory project using measured forces on hori- 
zontal submerged cylinders, in an infinite fluid, correlated 
with measured kinematics, Keulegan et al. (7) found that C, 
and Cr depended on a period parameter defined as =n, where 
Un is the maximum horizontal velocity. This parameter 
emphasizes the relationship of the force coefficients to 


the amplitude of the particle displacement and is related 
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to the ratio of the displacement to cylinder diameter ac- 


cording to Airy wave theory by 


ip 
A/D = 57 7 - (5) 


They found that with a given cylinder diameter, Cy varies 


as the maximum current, Uey is varied. A similar variation 


occurred for Ch: A correlation between these coefficients 


and the Reynolds number as found by Dean et al. (3) did not 
appear to exist. In this regard, Thirriot, Longree, and 


Barthet (13) found that when A/D > 10 the C, for steady 


D 
flow was approximately equal to the Ch for oscillatory flows 


for a vertical cylinder. Also for 1 S$ A/N £ 10, the Cy FOr 


oscillatory flow was larger than the C,, for steady flow. 


D 
This variation was valid for Re s 4 x ior. For larger Rey- 


nolds numbers the value of Cp for oscillatorv flow equalled 


the Ch for steady flow regardless of the particle displace- 


ment. 


Specifically, Keulegan et al. (7) showed that Cy = 1).9 


UmT UmT 
for = < 5 (A/D < 0.8) and increased to Cy = 2.5 at a eee 


15 (A/D = 2.4). For greater values of the period parameter 
there was a gradual decrease in the drag coefficient to the 


steady state value found by Roshko (9). They found that 


UmT 
‘a 2.0 fOr = < 6 (A/D < 1.0), as expected from potential 


U 
flow theory. For 


C 


T 
5 Spaacmel YN ial urge 631 IL Gy Baar 


Yamamoto, Nath, and Slotta (15) investigated the de- 





pendence of both C. and Ch on the relative distance of a 


i 






au 


horizontal cylinder from a plane boundary. For potential 
flow, where A/D < 1.0, they found that Cy = 4,15 at e/D = 
0.02 compared with Cr = 2.0 for a cylinder in an infinite 


fluid. They also reported that C, increases to large nega- 


L 
tive values as the cylinder is placed nearer a plane boun- 
dary for potential flow. They showed Ch, = -8.,45 for e/D = 
0.02. In addition, they emphasized the need for further ex- 
perimentation to determine the effect of variations in water 
depth on force coefficients, although they could not deter- 
Mine the effect because of limited experimental data. 

The experimental results discussed ahove show that the 
force coefficients for a circular cylinder do indeed depend 
on at least four factors: the Reynolds number, the 
Keulegqan-Carpenter period parameter, the proximity of a 
plane boundary, and variations in water depth. Values of 


these coefficients must therefore be chosen carefully when 


used in the Morison equation. 


Limitations of the Morison Fquation 


The Morison equation is limited in application and does 
not provide realistic values of forces in the following 
cases. First, the forces predicted by this equation are 
necessarily dependent on a theoretical wave theory if flow 
kinematics cannot be actually measured. In this regard some 
of the differences between two widely used wave theories 


will be explored in the next section. Second, according to 
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Ippen (6), the force on a cylinder in a wave pressure field 
as found by the Morison equation will not take into account 
the mean pressure surrounding the cylinder. Thus, the cyl- 
inder must be additionally designed to withstand the ex- 
pected hydrostatic pressure. Third, the assumption of con- 
Stant values of drag and inertia coefficients throughout the 
water depth and wavelength is consistent for cases where the 
flow kinematics are actually measured. Here the force co- 
efficients are dependent on Keulegan-Carpenter parameter. 
However for cases where theoretical kinematics are used, the 
drag coefficient varies with the Reynolds number while the 
inertia coefficient remains constant. Fourth, as discussed 
by Yamamoto et al. (15), the proximity of a plane boundary 
to the cylinder alters these force coefficients greatly. 
For example, according to potential flow theory, the small 
gap between a cylinder and the bottom boundary can produce 
large vertical forces toward the boundary, an effect which 
decreases as the cylinder is moved from the boundary. lLike- 
wise, the coefficients are affected by the proximity of a 
free surface. They also indicated that there may exist a 
phase shift between the ambient velocity and the drag force. 
Thus, the coefficient of inertia evaluated under the assump- 
tion of a zero drag force may be in error due to this phase 
Shift possibility. Fifth, for a cylinder resting on the 
bottom, the Morison equation predicts no vertical forces; 


whereas due to the asymmetry of the flow, a net positive 






Ihe 
vertical force occurs. Sixth, the force coefficients may be 
influenced by the water depth. And finally seventh, the 
Morison equation predicts no transverse forces due to vortex 


shedding or periodic wake phenomena. 


Comparison Between Airy Wave Theory and Stream 
Function Wave Theory 


In order to use an equation such as the Morison equa- 
tion for deSign purposes, values for the force coefficients 
and the fluid kinematics surrounding the structure must be 
obtained. The force coefficients are obtained experimental- 
ly either using a laboratory model or a field prototype. 

The experimentally measured forces are correlated with 
either measured or theoretical values of the fluid veloci- 
ties and accelerations. However, in general practice, fluid 
kinematics surrounding the structure are not known a priori. 
Generally only surface fluctuations have been recorded. 
Therefore, the use of a wave theory which will predict 
accurate values of the fluid kinematics is necessary to 
determine satisfactory values of the force coefficients. 

Although there are numerous wave theories available to 
the analyst, a comparison will be made between the first 
order Airy wave theory and the Stream function wave theory 
established by Dean (2). This brief comparison will empha- 
Size the differences in predicted surface profiles and kine- 


matics for both theories in order that these differences may 
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later explain some of the scatter found in force coefficient 
determinations in this investigation. 

In review, the Airy wave theory, used throughout this 
investigation to determine fluid kinematics, is the classi- 
cal hydrodynamic solution to the irrotational and incompres- 
sible progressive wave. It is the first order solution to 
the Laplace equation with boundary conditions evaluated at 
the still water level, where the Laplace equation for two 
dimensional flow is 


2) 
ay 


= 
According to the perturbation method of solution, higher 
order solutions to the Laplace equation will more nearly 
satisfy the free surface boundary condition. The Stream 
function theory is applicable to waves of steepness up to 
breaking. It also satisfies the two dimensional Laplace 
equation and assumes the wave flow does not change with res- 
pect to time for cases based on measured wave profiles. It 
is formulated by a nonlinear numerical perturbation proce- 
dure which reduces the error in the dynamic free surface 
boundary condition to acceptable levels. Since it satisfies 
the kinematic free surface boundary condition exactly due to 
the steady state assumption, it provides a sianificantly 
better fit to these boundary conditions for shallow and deep 
water waves than most analytical theories. Currently the 


Stream function theory has been evaluated for 40 different 
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waves varying from shallow to deep water conditions and from 


is the breaking wave 


B/E, = 0:25 to H/H, = 1.00, where H, 


height. 

In representing surface profiles of various waves, the 
Stream function generally predicts higher peak wave ampli- 
tudes and a more shallow trough than the symmetric predic- 
tion by Airy theory. Figure 3 shows the dimensionless sur- 
face profile according to the Stream function plotted aqain- 
St wave phase for the case where H/H, = 0.75, H/T? = 0,118 
ft/sec’, and h/t? = 0.20 ft/sec’. Figure 4 again shows the 
dimensionless surface profile predicted by the Stream fune~ 
tion compared to the experimental run 723 and the profile 
predicted by Airy theory. From these figures one would ex- 
pect that some of the actual kinematics of flow during the 
passage of a wave would be poorly predicted by Airy wave 
theory. However, a notable exception to this was pointed 
out by Grace and Rocheleau (5) and Dean and LeMéhauté (4). 
They indicate that Airy wave theory predicts the horizontal 
velocities of the flow exceptionally well near the bottom 
boundary beneath the wave crest. 

Figure 5 compares values of dimensionless horizontal 
velocity and acceleration and vertical acceleration for both 
wave theories. Since the cylinder is far from the free sur- 
face, the presence of the cylinder does not affect the flow 
field. Therefore these particle kinematics represent those 


which the cylinder experiences at @/D = 0.042. The maximum 






n/H 


Brgure 3. 
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Figure 5. Kinematic comparison for Airy theory and Stream 


function theory for cylinder at e@/D = 0.042. 
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value of the horizontal velocity and acceleration are nearly 
equal in both theories. In fact, for this example, the 
maximum horizontal velocity according to the Stream function 
is 16 percent greater than predicted by Airy theory. Like- 
wise, the horizontal acceleration according to the Stream 
function is only eight percent greater than predicted by 
Airy theory. A portion of this difference in acceleration 
is justified since the Stream function theory evaluates the 
total acceleration while Airy theory evaluates only the 
temporal acceleration. However, the horizontal acceleration 
predicted by the Stream function shows a definite phase 
shift toward the crest of the wave. The dimensionless 
vertical acceleration predicted by the Stream function ex- 
hibits a large difference from that predicted hy Airy wave 
theory at most wave phases. In fact, at the crest, the 
vertical acceleration according to the Stream function is 
over 160 percent greater than predicted values of Airy 
theory! Therefore, it may be concluded that maximum values 
of horizontal velocity and acceleration determined by Airy 
theory are fairly accurate under the crest of a wave and 
near the bottom. However, values of vertical acceleration 
may be seriously in error throughout the water depth for 


nonlinear waves. 
Conclusion 


The Morison equation not only fails in the prediction 


of certain forces on cylinders placed within one cylinder 
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diameter of a plane boundary but when the flow kinematics 
are theoretically determined, the force coefficients within 
that equation may vary throughout the wave. These coeffi- 
Cients depend on the Keulegan-Carpenter period parameter for 
measured kinematics and on the Reynolds number for theoreti- 
cal kinematics provided by the Stream function. Therefore, 
it was the general purpose of the research reported herein 
to investigate the effects of a plane boundary and water 
depth on the forces on a horizontal circular cylinder in 
oscillatory flow for as wide a range of Reynolds numbers and 
values of A/D as possible. 

Utilizing several results from Yamamoto et al. (15), a 
model is presented here which includes the basic approach of 
the Morison equation but also includes additionai factors 
which must be considered in the evaluation of forces on 
horizontal cylinders. This model will be used later to sum- 
marize some experimental results. For uniform flow, the 


horizontal and vertical forces can be considered as 


dU dU 
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where Cr = Inertia coefficient. For a uniform flow, 


C, - (1.0 + Cy) « 


Cu is an added mass coefficient which for a 
cylinder uniformly accelerated in a fluid at 
rest, is equal to the ratio of the hydrodynamic 
force necessary to maintain the motion of the 
cylinder in the presence of the fluid to the 
displaced mass of the fluid by the cylinder. In 
anwdnocinite fluid, Cy = 1.0 according to poten- 
tial flow theory. For the case of a stationary 
cylinder in a uniformly accelerating flow, a 
coefficient due to the pressure field causing 
acceleration of the fluid is added to the added 
mass coefficient to yield the inertia coeffi- 
cient. This pressure coefficient is always 
equal to 1.0 for any submerged object in a uni- 
form flow. It is noted here that for a uniform 
flow, variations in e/D will produce variations 


in Ga by varying C For a nonuniform flow the 


MW M° 


pressure coefficient no longer equals 1.0, nor 


does Cu remain defined as before. Thus, vari- 


ations in Cc, result from variations in e/D and 
flow nonuniformity. 
Cr = Inertia coefficient for convective acceleration 


which is defined similarly to the inertia 
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coefficient for local acceleration. From poten- 


tial flow theory, C, = 2.0 and this value is 


A 
constant for all values of e/M. It is noted 
here that values of Cy determined by Dean et al. 
(3) using the Stream function wave theory were 
based on total fluid acceleration including the 
convective contribution. 


£ (6) 


Phase shift factor between ambient flow and the 
drag force which is to be experimentally 
determined. 

[T = Induced circulation around the cylinder due to 
vortex formation which is combined with the uni- 
form velocity in the well-known Kutta-Joukowski 


relation. 


The horizontal convective acceleration of the ambient 


Flow 


is constant throughout the water depth according to Airy 
Wave theory and varies at twice the wave frequency. The 


vertical convective acceleration of the ambient flow 


oV oV 

OX Ch 
is non-periodic and always in the positive direction. For 
most of the experimentation reported herein, the horizontal 


and vertical convective acceleration as calculated by Airy 
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wave theory was less than ten percent of the temporal ac- 
celeration in those directions. The convective contribution 
to the total acceleration was consequently ignored in the 
determination of the inertia coefficient. 

Therefore, with the conclusion of this review, the 
mechanics of unsteady flow across a cylinder and the effect 
of a plane boundary on the hydrodynamic forces on the cylin- 


der will be examined. 
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IIIT. EXPERIMFNTAL INVESTIGATION 


Apparatus 


Wave Research Facility 


Experiments were conducted at the Oregon State Univer- 
Sity Wave Research Facility during August and September, 
1974 to further investigate the effect of a plane boundary 
on a horizontal circular cylinder placed in wave flow. As 
shown in Figure 6 the wave basin is 104.27 m (342') long, 
3.66 m (12') wide, and generally 4.57 m (15') deep and pro- 
vided sufficient length and depth possibilities to generate 
large amplitude waves with minimum beach reflection. Beach 
configuration was adjusted to a 1:12 slope which provided a 
maximum reflection coefficient of eleven percent and a mini- 
mum of three percent with an average value of seven percent. 
The wave generation system is hydraulically activated and 
electronically controlled from a station near the test cylin- 
der location. The system was capable of producing clean, 
uniform waves with wave heights up to 1.52 m (5'). 

Bottom configuration was established in Figure 6. Each 
3.66 m (12') square concrete slab was individually lifted 
to provide the needed false bottom and beach configuration. 
Fach was secured in place by heavy clip angles bolted to the 
Sides of the wave basin. The three slabs within the test 


area were carefully sealed to prevent uplift currents due to 
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the false bottom from altering the flow around the cylinder 
with the passage of a wave. The edges were sealed with 
Strips of treated 1/2" plywood while the spaces between 
slabs were sealed by 1/2" x 6" x 10' lengths of flat steel 
plate. The test holes were sealed with No. 9 rubber 


stoppers. See Figure 7. 


Test Cylinder 


The apparatus itself was designed to fulfill the re- 
quirements of the investigation. It was a cylinder capable 
of transmitting wave forces in terms of deflection or strain 
to recording instruments. It was to remain rigid during the 
passage of the wave, but was also adjustable to various 
heights above the bottom. 

As shown in Figures 7 and 8, the apparatus consisted of 
a three-piece 12" OD x 1/4" 6061-T65il1l extruded aluminum 
tube where the center test section was 1.07 m (3.5') long 
with each of the two end sections at 1.22 m (4') long. The 
two end sections, or dummy sections, were required to ap- 
proximate two-dimensional flow conditions around the test 
cylinder. The spacing between the test section and the end 
sections was from 3 mm to 5 mm. Two threaded 1" diameter 
Stainless steel split-type threaded shaft collars were 
placed above and below the dummy section at each steel bar 
to hold the cylinder and steel bar rigidly. Steel guy wires 


extended in the direction of flow from ahove and below each 
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29 
dummy section at each steel bar to points of attachment on 


the test Slab. 


Force Dynamometer 


In order to measure deflections of the test cylinder 
from pressure distributions at its surface, a strain gage 
dynamometer was designed. Two 1" diameter 6061-T6 aluminum 
bars were machined as in Figures 9 and 10 to accommodate 16 
Micro-Measurements CEA-13-189UW-120 strain gages. The gages 
measured horizontal and vertical deflections of the test 
section and were attached as in Figure 1l. Each end of the 
test section and the adjacent ends of the dummy sections 
contained a rigidly fastened box of 1/2" aluminum plate 
through which the strain gage bars passed and were rigidly 
fastened. Thus, at each end of the test section the strain 
gage bars were represented as a beam fixed against all 
motion at one end and constrained to only horizontal and 
vertical displacement at the end within the test cylinder. 
The test structure was designed to have negligible deflec- 
tion with respect to the diameter. All strain gages were 
water-proofed. 

Output from the strain gages were conditioned and 
balanced through four Techtronix 3C66 Carrier Amplifier 
Bridges. Two sonic profilers were positioned above the 
cylinder and forward of the cylinder (toward the wave board) 


to measure the water surface fluctuations. The profilers 
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bi! 
were 3.94 m (12.93') apart for runs 107 to 1430 and 4.44 m 
(14.58') apart for runs 1501 to 1617. Conditioned signals 
from the strain gages and signals from both profilers were 
recorded on a twelve-channel Honeywell 1508 Visicorder. The 
signal from the sonic profiler over the cylinder was also 
recorded at the wave generating station by a Hewlett-Packard 


two-channel chart recorder. 


Test Procedure 


A series of four water depths were deSigned in an array 
with six gaps between the cylinder and the plane boundary. 


Values for water depths and gaps were: 


emt) e [cm (in) ] 
Ondo (2.35) 0 

0.957 (3.14) Vay) COS) 
Oe (52) Selo (22 5) 
Zia (7.8) Slice (20) 


be 24 (620) 
SOR ae (12,0) 


60-96 Sa022. 0) 


Of these 24 combinations, 16 were chosen to he investigated. 

Each of the 16 combinations contained individual runs of 

varying wave periods and wave heights. Specific ranges were: 
te? 40, 96450 6.20 sec 
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34 
H/T2 : 0.0085 to 0.505 ft/sec2 


2 . 0.078 to 2.57 £t/sec2 


h/t 
Re : 2.0 x 103 to 2.0 x 10° (based on maximum hori- 
zontal velocity). 
Thus, the experimental schedule was designed to establish 
information on the effect of a plane boundary on wave forces 
on a horizontal cylinder whose longitudinal axis was perpen- 
dicular to the direction of flow. 

The test procedure covered the following steps: 

1. The test cylinder was lowered or raised to the 
position of the required experimental configuration. Move- 
ment of the cylinder was controlled by a chain hoist and 
Spreader bar for slow, even adjustment so as not to over- 
Strain the strain gages. With the proper configuration, 
shaft collars were secured and guy wires adjusted to main- 
tain a rigid structure. 

For the cylinder at e = 0, strips of foam rubber 2.5 cm 
(1") thick were placed between the test cylinder and the 
bottom as a flow barrier with the cylinder and the bottom 
as a flow barrier with the cylinder secured at the position 
of e = 1.27 (0.5"). For runs 6 and 12, the foam rubber bar- 
rier was also placed between the dummy sections and the bot- 
tom. The resistance of the foam rubber to vertical loading 
of the test cylinder was assumed negligible and no external 
Calibration of the cylinder was conducted for this confiaqu- 


ration. Therefore, the only difference between the set up 
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of configurations for e = 0 and e = 1.27 cm (0.5") was the 
insertion of the foam rubber barrier. The adjustments to 
the test cylinder were all conducted while submerged using 
two personnel in scuba gear. Each adjustment generally took 
30 minutes time. 

2. The wave generator was set to the required frequen- 
cy and amplitude. Check runs were required occasionally to 
assure the generation of the proper wave. 

3. With still water in the wave basin, each pair of 
strain gages were balanced and internally calibrated prior 
to each run. 

4. The test wave was then generated and results from 


the strain gages and sonic profilers were recorded. 


The above steps were repeated for each run at each re- 


quired water depth. 
Calibration Procedure 


The response of the test cylinder to horizontal and 
vertical forces was calibrated four different times during 
experimentation. As shown in Figure 12, known weights were 
attached to the test section via a minimum friction pulley 
System in the negative horizontal and positive vertical 
directions. Calibration of the test cylinder was not con- 
ducted in the positive horizontal and negative vertical 
directions because of the symmetry of the force dynamometer 


system. The microstrains recorded were then correlated to 





Figure 12. 
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Horizontal and vertical calibration setup. 
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oY 
the known forces. The correlation remained nearly constant 


throughout the testing schedule. 


Natural Frequency of system and Resonance 


At the conclusion of all scheduled test runs the natu- 
ral frequency of the structure was roughly determined by re- 
cording the frequency of vibration after giving the test cy- 
linder an initial impulse. The following results were 
obtained: 

1. For cylinder in air, guyed above and below, with 

Seah 9orcm (24 is Sf.3= 23 HZ and fe sl Ae id 9 92 


V 
where fy and fy are the vertical and horizontal 
natural frequencies of vibration, respectively. 
2. For cylinder in air, guyed above, with e = 15.24 cm 

USA E f. = 23 Hz and thy =o! Se Omiriz 
Ba For cylinder in water at depth h = 1.07 cm (3.5'), 

guyed above and below, with e = 60.96 cm (24"): 
fy = Gg ,00hz ang fu =a" 52. HZ. 

The natural frequency thus determined is sufficiently 
greater than expected frequencies of transverse vibrations 
due to vortex shedding such that a resonance condition is 
unlikely to occur. This can be shown with the crude use of 
the Strouhal number which is a similarity parameter ex- 
pressed as: 
£ D 
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where f. = frequency of transverse vibrations due to vortex 
shedding. If S = 0.21 over the range 10° < Re < 10” and 
U = 1.07 m/sec then f . = 0.7 Hz for this test cylinder. 


Thus, for the case where e/D = 2.0, 


= _ ia lg >> ne 






IV. RESULTS AND DISCUSSION 


- Calibration Factors 


Results of the calibration procedures described in 
Chapter III are plotted in Figures 13 and 14. A linear re- 
gression procedure was applied to obtain a best fit line be- 
tween the values of microstrain and calibration loading. 
This relationship yielded the calibration factor for both 
the horizontal and vertical forces such that the microstrain 
recorded from the movement of the test cylinder when multi- 
plied by this factor yielded the total load on the test 
cylinder. The resulting calibration factors used through- 


out the data reduction procedure were: 


OQ 
I 


= 0.252 pounds/microstrain 


QO 
; 


= 0.189 pounds/microstrain 


Wave Force Data Reduction 


Horizontal and vertical wave force traces and two wave 
profile traces were reduced for each run to yield the force 
data found in Appendix B. Wave celerity and observed wave 
length were calculated for most runs utilizing the known 
distance between the two sonic profilers. 

Raw data for each run were composed of six traces, four 
Of which were the variations of strain in the horizontal and 


vertical directions for both the left (west) and right 
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42 
(east) ends of the test cylinder. This duplication of 
strain readings provided a check on the alignment of the 
force rod at either end of the test cylinder. The forces 
were equal at both ends of the cylinder. However, the four 
gages on the right end of the test cylinder which recorded 
vertical deflections gave signals which continually drifted 
and prohibited determination of vertical forces. Thus, the 
results from the eight left end gages were used to provide 
the horizontal and vertical forces listed in Appendix B. 
The remaining two traces were the wave profile forward of 
the cylinder and above the cylinder. 

Appendix A contains the results of the computer routine 
which calculated various flow parameters for each run 
according to Airy wave theory utilizing experimental wave 
profile data. For all cases, the wave length calculated 
from Airy wave theory was utilized in lieu of observed 


values of wavelength. 


Flow Mechanics and Effect of Plane Boundary 
and Fluid Displacement on Forces 


As discussed in Chapter II, the hydrodynamic forces on 
a horizontal cylinder during the passage of a wave vary with 
the proximity of the plane boundary, expressed as e/D, and 
the relative horizontal displacement of the water particles, 
expressed as A/D. Utilizing experimental results and the 


analytical and empirical considerations of Chapter II, the 
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mechanics of flow around a horizontal cylinder will be des- 
cribed for three cases. 

Case 1 will consider the flow around a horizontal cylin- 
der sealed to the plane boundary. Runs 605, 612, and 620 
will be discussed for which e/D = 0. 

Case 2 will consider a horizontal cylinder near the 
plane boundary so that the definite effects of the plane 
boundary on the forces will be determined. Runs 705, 709, 
718, and 725 will be discussed for which e/D = 0.042. 

Case 3 will exhibit the flow around a cylinder in an 
infinite fluid. Essentially this will mean that the cylin- 
der will experience no effects from a plane boundary or free 
surface. As shown by Wilson et al. (14) and Yamamoto et al. 
mo), the lift force on a horizontal cylinder for e/D = 1.0 
is five percent to fifteen percent of the drag force with 


4 to 10°. Runs 1110, 1118, and 1120 


Reynolds numbers of 10 
where e/D = 1.0, will be utilized to describe this case. 
Thus, the location of the plane boundary and the tran- 
Sition from potential flow to the characteristics of real 
flow produces important variations in the types and direc- 
tions of forces on the horizontal cylinder. As found in 
Chapter II, one of the most useful parameters for locating 
transition points is A/D, the ratio of the maximum horizon- 
tal particle displacement to cylinder diameter. The value 
of the maximum horizontal particle displacement utilized 
throughout this investigation was calculated by Airy wave 


theory. 
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Case 1 - Cylinder Resting on Plane Boundary where e/D = 0 


For a cylinder resting on and sealed to a plane boun- 
dary, it is well known that the horizontal velocity will re- 
sult in horizontal forces due to the pressure drag of the 
wake when the wake develops and positive verticai forces due 
to the asymmetry of the flow around the cylinder. Also, the 
horizontal acceleration will produce horizontal inertia 
forces. The vertical velocity and acceleration are very 
small near the bottom and have a negligibly small effect on 
the cylinder compared with other forces. Figures 15, 16, 
and 17 show the experimental results for runs 605, 612, and 
620 where values of A/D are 0.09, 0.51, and 3.20 
respectively. 

From examination of all data for e/D = 0, the following 
general tendencies were found. For A/D < 0.3, the data show 
the horizontal and vertical force variations remain nearly 
Sinusoidal and increase linearly with wave height. Thus, at 
small values of particle displacement the flow is essential- 
ly potential flow and both horizontal and vertical forces 
are due to particle acceleration, as exemplified by the re- 
Seecs for run 605 in Figure 15. For A/D > 0.3, as shown in 
Figure 16, the vertical force variation ceases to be sinu- 
Soidal and begins to deform to exhibit two positive and two 
negative force peaks. It is at this value of A/D that 
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force peaks increase proportionally to the wave height 
squared and this parabolic increase verifies that they 
represent a poSitive lift force generated by the asymmetri- 
cal flow across the cylinder. The horizontal forces remain 
nearly sinusoidal in variation for all values of A/D investi- 
gated with deviations attributed to wave nonlinearity. The 
maximum horizontal forces remained larger than the maximum 
vertical forces. 

Figure 18 represents a summary of the foregoing discus- 
Sion. It shows the temporal variation of the wave surface 
with the expected positions of the maximum values of horizon- 
tal and vertical velocities and accelerations according to 
Airy wave theory. The relative magnitude of the experimen- 
tal values of the horizontal and vertical forces of the 
cylinder at various wave phases are also shown schematical- 
ly. From this the inertia dominated flow at low values of 
A/D can be seen as well as the onset of the horizontal drag 
force due to wake effects at large A/D values. The negative 
vertical force under the node of the wave increases linear- 
ly with wave height and thus is proportional to the vertical 
acceleration. The location of the force at a wave phase 
where, according to Airy wave theory, the vertical accelera- 
tion is zero, exhibits the inaccurate prediction of the 
vertical acceleration by that theory. As was shown in 
Figure 5, the difference between v according to Airy theory 
and v according to the Stream function is quite marked in 


both magnitude and direction for nonlinear waves. 
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Case 2 - Cylinder Near Plane Boundary where 0 < e/D < 1.0 


Tt has been shown that the lift force due to flow asym- 
metry is always upward for both potential flow and wake flow 
when e/D = 0. The transition from potential flow to wake 
flow is less clear for the case where a cylinder is near a 
plane boundary. As mentioned in Chapter II, one can now ex- 
pect negative lift forces to occur from the horizontal ve- 
locities due to potential flow effects. It will be shown 
that positive lift forces also occur due to horizontal wake 
effects. The remaining inertia and drag forces in the 
horizontal and vertical directions occur as the result of 
corresponding accelerations and velocities similar to Case 
1. Figures 19, 20, 21, and 22 give experimental results for 
runs 705, 709, 718, and 725 where values of A/D are 0.09, 
0.41, 1.62, and 3.43 respectively. 

From examination of all data for which 0.04 S e/D & 
0.5, the following observations were made. A small particle 
displacement again indicated that a wake has little chance 
of developing around the cylinder and consequently the 
forces on the cylinder can be discussed according to poten- 
tial flow considerations. For A/D < 0.25, the horizontal 
and vertical forces vary nearly sinusoidally and result from 
an inertia dominated potential flow as in Figure 19. For 
larger values of A/D where 0.5 < A/D < 1.0, and still with- 
in potential flow as in Figure 20, the horizontal velocity 


contributes to negative vertical forces since the higher 
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Figure 19. Horizontal and vertical forces, e/N = 0.042. 
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forces, e/D = 0.042. 
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Figure 21. Horizontal and vertical forces, e/D = 0.042. 
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velocity flow between the cylinder and the boundary gener- 
ates a lower pressure in this area. These observations can 
be verified. For A/D < 0.25, -F, occurring near the wave 
crest is proportional to wave height and for greater parti- 
Cle displacement, the variation in “Fy becomes nonlinear 
increasing with respect to H and proportional to me 

On the other hand, a large particle displacement im- 
plies that a wake will form and the forces on a cylinder 
must now include real flow or wake effects. In Figures 21 
and 22, the horizontal forces become a combination of in- 
ertia and drag forces but less obvious changes occur with 
the vertical forces. In Figure 21 for which A/D = 1.62 
the negative vertical force peak under the crest of the wave 
has shifted forward of the crest and a positive vertical 
force exists under the trough of the wave. Thus, the posi- 
tive vertical force is a lift force due to the wake behind 
the cylinder as in steady flows. The existence of this 
positive lift force is explained as follows. The wake slows 
the flow between the cylinder and the boundary from the in- 
crease in negative pressure behind the cylinder resulting 
from the clockwise rotation of the lower point of flow 
separation on the cylinder. Higher flow is readjusted over 
the cylinder and the pressure gradient becomes negative up- 
ward. This wake effect may be referred to as a "hydraulic 


Darrver,." 
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The concept of a "hydraulic barrier" is substantiated 
in Figure 23. This figure compares the positive lift forces 
discussed previously for e/D = 0 to the positive lift forces 
for e/D = 0.042. As explained in Chapter III, the change 
from the experimental configuration for the 700 numbered 
runs to the 600 numbered runs was made with the insertion of 
a foam rubber barrier between the cylinder and the bottom. 
The positive vertical force increases parabolically with H 
for both cases of e/D = 0 and e/D = 0.042. The similarity 
in the positive vertical forces suggests that the "hydraulic 
barrier" acts in a manner similar to the foam rubber bar- 
rier. The negative vertical force for e/D = 0.042 first in- 
creases parabolically with H and approaches a constant value 
for A/D > 1.5, whereas the negative vertical force for e/D = 
0 stays negligibly small for all H. Therefore, for the case 
where 0 < e/D < 1.0 and A/D > 1.5, the cylinder experiences 
large alternating lift forces occurring at twice the wave 
frequency. This suggests a possible failure mode for sub- 
marine pipelines. 

Figure 24 presents a summary of the foregoing discus- 
sion in a format similar to Case 1. Again, as discussed in 
Case 1, the negative vertical force under the node of the 
wave, in contradiction to Airy wave theory, may be due to 


i?) 
~V if a higher order wave theory is used. 
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Case 3 - Cylinder Far from Plane Boundary and Free Surface 
where e/D 1-20 


Figures 25, 26 and 27 for runs 1110, 1118, and 1120 are 
presented for this case where the cylinder is far from the 
bottom boundary and free surface. Values of A/D are 0.75, 
1.72 and 3.18, respectively. 

From examination of all data for which e/D 2 1.0, the 
following tendencies were noted. For A/D <* 2.0, the nega- 
tive vertical force peak near the wave crest increases in 
proportion to H. For 2.0 £ A/D § 3.0, -F,, increases non- 
linearly with H and at A/D > 3.0 this rate of increase slows. 
This trend for vertical forces is similar to that of Case 2, 
but the transition occurs at different A/D values. At 
A/D = 3.0 the wake should already be well developed such 
that the decrease in -F,, cannot be completely explained by 
assuming a transition from potential to real flow, unless 
values of A/D calculated according to Airy theory are in 
error. In addition, there exists no appreciable phase shift 
of the -F,, peak forward of the wave crest as was observed in 
Case 2. This observation is expected since there exists no 
appreciable positive lift force in the absence of a plane 
boundary. In Case 2, the generation of the positive lift 
force at the wave crest resulted in shifting the negative 
vertical force peak forward of the wave crest. 


The vertical forces vary with one negative and one 


positive force peak for A/D < 0.8 characteristic of an 
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Figure 26. Horizontal and vertical forces, e/D = 1.0. 
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Figure 27. Horizontal and vertical forces, e/D = 1.0. 
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inertia dominated potential flow. For larger values of A/D 
up to the limit of this investigation at A/D = 5.0, the 
vertical forces are characterized by a two- or three-cycle 
alternating negative and positive force peak within each 
wave length. 

Figure 28 summarizes this discussion. For run 1120 for 
instance, the alternating positive and negative vertical 
forces are not adequately predicted by the Airy theory kine- 
matics. Use of the Stream function for prediction of verti- 
cal accelerations may explain some discrepancies similar to 
the previous cases. As shown in Figure 29, the variation 
in the vertical force for Run 1125 shows a close correlation 
to the variation of the vertical acceleration predicted by 
the Stream function wave theory. However, there appears to 
be a need for further investigation into the problem of 
periodic vertical forces which are not explained by either 
inertia or drag effects. It is only suggested here that 
these forces may be the result of the superposition of in- 
ertia and drag effects and the alternating vertical pressure 


effects due to vortex shedding. 


Summary and Conclusion 


In summary, one can see from the discussion of the 
three cases that the maximum forces experienced by the 
cylinder depend on the proximity of the plane boundary and 


the kinematics of the flow field. Table 1 is presented as 
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Figure 29. Vertical force on horizontal cylinder correlated 
with vertical acceleration computed by Stream 
function theory. 
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a summary of expected effects from potential flow and real 
flow depending on the location of the cylinder with respect 
to a plane boundary. The table gives the contribution of 
the horizontal and vertical forces to the total force on the 
cylinder in terms of Equations 6 and 7. Contributions due 
to convective acceleration were included when the convective 
acceleration was at least ten percent of the local accelera- 
tion as calculated by Airy wave theory. It also lists the 
location of the maximum horizontal and vertical forces with 
the phase angle measured positive in the direction of the 
wave advance with ? = 0 at the wave crest. Finally, the ex- 


pected range of the ratio of maximum forces, F is entered 


MM’ 
such that 


F 

Hmax 
F.. ==>. 
M F 

Vmax 


Evaluation of Force Coefficients 


In order to evaluate the various force coefficients of 
inertia, lift, and drag as found in Equations 6 and 7, 
values of pertinent wave forces were compared with corre- 
sponding velocities and accelerations predicted by Airy wave 
theory. In each case a linear relation was expected and 
the slope of the line which best visually fit the data 


yielded the coefficients. Results are given in Table 2. 
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Inertia Coefficient 


For small values of A/D, horizontal and vertical forces 
on the cylinder will generally be due to inertia effects and 
proportional to U Or V if convective acceleration effects 
are ignored. Figures Cl to C9 compare to ae near the wave 
crest with V evaluated at the force peak where A/D < 1.5. 
The linear relationship between the inertia force and ac- 
celeration determine the value of Cy as defined in Equation 
2. Nonlinearity and scatter occur for at least two reasons. 
First, the effect of the plane boundary increases “Fy para- 
bolically with respect to H for potential flow. And second, 
the value of the vertical acceleration calculated by Airy 
wave theory does not include convective acceleration ef- 
fects. The use of the Stream function theory which predicts 
values of the total acceleration may reduce some of the 
scatter. See Chapter II. 

In Figures C10 to Cl7, TP naxx located between the node 
behind the wave crest and the trough is plotted against v 
evaluated at the force peak for A/D < 1.5. The values of 
Cr at this position of the wave were generally lower than 
those values at the wave crest. 


In Figures C18 to C28, -F located. near the node be- 


Hmax 
\) 
hind the wave crest is plotted against URES. The maximum 
o 
value of U was utilized throughout since according to 


Figure 5, a comparison between Airy and Stream function 
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theories shows maximum values are nearly equal for both 
theories. 

In Figures C29 to C36, +F located near the node 


Hmax 


co) 
forward of the wave crest is plotted against U It is 


max- 
apparent that nonlinearities occurring here may be due to 


the residual drag forces which increase +F,, parabolically 


H 
with respect to H. Again, values for C. evaluated from 
both the vertical and horizontal forces are plotted against 
values of e/D in Figures 30 and 3l respectively. From 
Equation 8 the variation in the coefficient of inertia for 
a Circular cylinder is essentially the variation in the 


added mass coefficient. This verifies experimental values 


found by Yamamoto et al. (15) and Schiller (12). 
Lift Coefficient 


It has been shown that one of the effects of a plane 
boundary on the forces on a horizontal cylinder is to create 
both positive and negative vertical forces which are propor- 
tional to the horizontal velocity squared and thus called 
lift forces. With these forces separated from the total 
force and plotted against vu, the positive and negative 
values of C,, can be determined from Fquation 4. 

To separate the negative lift force from the total 
negative vertical force at the crest of the wave, the as- 


Sumption that C. is constant throughout most of the wave 


uy 


is utilized. Since values of ie are generally much greater 
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Figure 30. Vertical inertia coefficient for circular 
cylinder near a plane boundary. 
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Figure 31. Horizontal inertia coefficient for circular 
cylinder near a plane boundary. 
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72 
than values of vm slight variations in values of Cr will 


negligibly affect values of C Using the maximum value of 


Lice 
Cy applicable to the run of interest, the negative lift 
force equals the total negative vertical force minus the 


negative vertical acceleration force, or 


where v is evaluated at the force peak. Thus separated, 
the negative lift force is plotted against te evaluated 
at the wave crest as in Figures C37 to C43. Deviations from 
the linear relation occur most frequently for runs where 
A/D > 1.5 and are apparently due to the increased wake ef- 
fect and the consequent generation of positive lift accord- 
ing to the "hydraulic barrier" concept developed in Chapter 
ny, 

The positive lift forces due to flow asymmetries can 


be readily located. For runs where e/D = 0, the +Fy _, 


2 and 


near the wave crest is consistently proportional to H 
can then be plotted against u- evaluated at the force peak 


to give values for C, as in Figures C44 to C51. A plot of 


L 
+E max at the trough of the wave against i evaluated by 
Airy wave theory exhibited more scatter since the horizontal 
velocity at the trough is generally less than at the crest 
for nonlinear waves. For runs where e/D > 0, the Eas 
which begins to appear near the wave crest as A/D increases 


is plotted against u- evaluated at the force peak. 
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Values of C. for both potential flow and real flow are 
plotted against e/D in Figure 32. Experimental values from 


Yamamoto et al. (15) and Schiller (12) are also shown. 


Drag Coefficient and Phase Difference Between Drag Force 
and Ambient Flow 


AS indicated by Yamamoto et al. (15), there exists a 
definite phase difference between the ambient fluid velocity 
and the drag force. For laminar flow they showed that the 
total drag force on a harmonically oscillating circular 
cylinder is 1/4 out of phase with the cylinder velocity. 
Therefore there remains a residual drag force on the cylin- 
der when the cylinder reaches maximum acceleration. For 
turbulent flow, this phase lag will also exist but at a 
value still to be determined. To illustrate that this 
phase difference in fact may exist for unsteady flow, 


Figure 33 compares values of -F,, and wave height for runs 


H 
1611-1615. The values of the horizontal forces at various 
phase angles from the crest show that a linear relationship 
occurs at 7/6 radians from the crest. Previous evaluations 
ona Cr assumed that a linear relation between force and wave 
height occurred at the node of the wave for all values of 
A/D; however the linear relationship at the node is ter- 
Minated for A/D > 1.5. For A/D > 1.5 linearity occurs at 
n/3 radians from the node. Thus, the assumption of zero 


drag forces at the node in determining Cr may be in error 


resulting in the underevaluation of Cy. 
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Figure 32. Lift coefficient for circular cylinder near a 
plane boundary at A/D < 1.5 and A/D > 1.5. 
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a7 
The value of the positive horizontal force occurring 


under the wave crest was utilized in evaluating C With 


1 


this force assumed to be F the coefficient of drag was 


ID. 
evaluated according to 

F 

C= D 


2 lk oDUS 


Max 


against oe ‘ 
ma 


Meoures C52 to C56 plot F - 


D 


In Figure 34 values of C. are plotted against e/D at 


D 
an average value of Reynolds number of 1.3 x aie Although 
considerable experimental investigation remains to verify 


this relationship between C. and e/D for various values of 


D 
Reynolds numbers, a comment about the peak value of C, at 


D 
e/D = 0.17 is necessary. 

In a current investigation by Hafen (17), using flow 
visualization techniques for steady flow across a horizon- 
tal cylinder near a plane boundary, it was shown that the 
width of the wake behind the cylinder, expressed as the 
distance between free streamlines, decreased with increas- 
ing values of e/D from 0 to 0.9 for a constant subcritical 
Reynolds number. For e/D > 0.9, the width of the wake in- 
creased with increasing values of e/D to the stable wake 
Size in an infinite fluid. The wake width at e/D = 0.9 was 
less than the wake width in an infinite fluid. In addition, 
vortex shedding became visible for e/D > 0.25 and subcriti- 
Cal Reynolds numbers. For e/D < 0.25, vortex shedding was 


Suppressed and the wake appeared uniform behind the cylinder. 
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In comparing the width of a wake for a circular cylin- 
der at subcritical Reynolds numbers to the drag coefficient, 


Roshko (16) found that for a given cylinder C. was inversely 


D 
proportional to the wake width. Thus, for subcritical Rey- 
nolds numbers one would expect an increase in Ch For a de- 
crease in wake width. As with the results from Hafen (17) 
one would expect an increaSe in Cy for increasing e/D values 
oe to 0.9. 

Considering the characteristics of steady flow past a 
cylinder near a plane boundary found above, it is possible 
that the increase in Ch for a cylinder near a plane boundary 
in unsteady flow may be attributed to the decrease in wake 
width. In addition, it is noted that the value of Cy for 
unsteady flow occurs near a value of e/D where vortex 
shedding begins to appear in steady flow. However, since 
no evidence of vortex shedding was found, it is difficult 
to correlate the mechanism of vortex shedding to the termi- 
nation in the increase of Cp: Additional experimentation 
and research is definitely required. 

Therefore the effect of the plane boundary on the drag 
coefficient for a horizontal circular cylinder in unsteady 
flow is that CD nearly tripled the value of CD in an in- 
finite fluid as the cylinder approached the plane boundary. 
Also for the cylinder sealed to the plane boundary, the Cy 


was 30 percent greater than the Cy in an infinite fluid. 
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BETCCEMOSENacer MeCEnson FOrce Cocrrtci1encs 


The evaluation of force coefficients utilized in re- 
lationships such as the Morison equation depends upon the 
assumption of a uniform flow field at the structure. If the 
size of the structure is large compared to a characteristic 
dimension of the flow, e.g. the wavelength or water depth, 
then the assumption of uniform flow is in serious error. 
merea horizontal circular cylinder in a fixed position with 
respect to a plane boundary, the effect of variations in 
water depth is to create a nonuniform flow. Thus variations 
in water depth produced variations in the coefficients of 
inertia and lift for A/D < 1.5. Water depth had no apparent 
effect on the coefficients of lift and drag for A/D > 1.5. 

Specifically, referring to Table 2, values of water 
depth with respect to cylinder diameter, expressed as h/D, 
femme. 35, 3.14, 5.32, and 7.8 for @€/D = 0.042 (runs 2, 5, 

7 and 13, respectively). Theoretically, for e/D = 0.042 and 
a uniformly accelerated flow, C. = 2.8. Values of C, from 


ih a 


vertical forces for this value of e/D varied from Cy = 1.82 


to Cy = 3.08 as h/D increased from 2.35 to 7.8. - Values of 


Cr from horizontal forces varied from Cr =—e2 5c 2cO Cr = 5) (UNS 


as h/D increased from 2.35 to 7.8. Theoretically, for e/D = 


0.042 and a uniform flow, C, = -7.0. As h/D increased from 


L 


ee 2.35 to h/D = 7.8, for A/D < 1.5, C varied from 


L 
Cc. = -4,29 to Cy, = eae ZO Therefore as water depth becomes 
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more shallow, the assumption of uniform flow becomes less 


valid for evaluation of both vertical and horizontal Cy and 


C,. where A/D < 1.5. However, the horizontal Cr is less af- 


fected by changes in water depth than the vertical Cr: Ee 
is noted that for h/D = 7.8, values of both vertical and 


horizontal C. are approximately equal and close to predicted 


r 
potential flow values. Consequently at h/D = 7.8 the effect 


of water depth on force coefficients appears negligible. 


In conclusion, the values of C. and Cy, as predicted by 


iE 
potential flow theory are conservative when applied in 
shallow water depths were h/D < 7.8. Also, variations in 
water depth affected vertical inertia coefficients more than 
horizontal inertia coefficients. Variations in water depth 


did not appreciably affect values of Cc, and CC. for A/D > 


D 
> 
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V. SUMMARY AND CONCLUSIONS 


From the experimental results, the following summary 
is made and conclusions drawn: 

1. The hydrodynamic force coefficients of inertia, 
lift and drag for a horizontal circular cylinder varied 
according to the proximity of the plane boundary, expressed 
as e/D, the relative horizontal particle displacement, ex- 
pressed as A/D, and the relative water depth, expressed as 


la D.. 


Effects of Particle Displacement 


2. For 0 < A/D < 1.5 the flow around the cylinder was 
that of potential flow to which potential flow theory is 
applicable. For A/D > 1.5, wake effects become significant 
and theoretical results from potential flow theory need to 


be modified. 
Effects of Plane Boundary 


Sm Values of C. from horizontal forces increased from 


C, = 36 to Cr = 3.47 as the cylinder approached the plane 


boundary. 
4. Values of Cy from vertical forces increased from 
C. = 1.92 to C, = 2.37 as the cylinder approached the plane 


1 2 
boundary. 
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5. The inertia coefficients obtained from the horizon- 
tal force data were nearly equal, but generally greater than 
those obtained from vertical forces. The experimental 


values of Cy agreed well with theoretical values for a uni- 


formly accelerated flow except for shallow water depths as 
noted below. 


6. For A/D < 1.5 and e/D > 0, values of Ch varied from 


Cy. = 0 to Ch = -5.3 as the cylinder approached the plane 


boundary. Experimental values agreed well with theoretical 
values for uniform flow except for shallow water deptns as 
noted below. 


g. For 1.5 < A/D < 5.0 and -é/D = 0, Ch = 190% 


MEE Gr), 59< A/D <. 5. 7/7and e/D > 0, Cy, increased from 


C., = 0 to CL. = 2.6 as the cylinder approached the plane 


boundary. 

9. The horizontal cylinder experienced alternating 
positive and negative vertical forces, which at times were 
of the same magnitude as horizontal forces, when 1.5 < A/D < 
5.7 and 0 < e/D < 1.0. These lift forces occurred at twice 
the wave frequency and may constitute a possible failure 
mode for submarine pipelines. 

10. For an average Reynolds number of 1.3 x LO? as the 
cylinder approached the plane boundary, Ch increased from 


Cy = 0.66 at e/D = 1.0 to Ch = 1.9 at e/D = 0.17, and then 


decreased to Cy = 0.93 at e/D = 0. 
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Effect of Water Depth 


ll. Relative variations in water depth, expressed as 


h/D, produced variations in C. and Ch: 


A 


12. As h/D varied from h/D = 2.35 to h/D = 7.8, Cr 
evaluated from vertical forces at e/D = 0.042 varied from 


m= 1.8 to C. = 3.08. 


i a 
13. As h/D varied from h/D = 2.35 to h/D = 7.8, Cy 

evaluated from horizontal forces at e/D = 0.042 varied from 
Cc. = ye, tO C,= 3.05. Variations in the horizontal force 
coefficient were less than the vertical. 

mae For A/D < 1.5 and e/D = 0.042, as h/D varied from 
[mee = 2.35.to h/D = 7.8, Cc, varied from Ch = -4,.29 to Cy, = 
=5 (Ae 

15. For A/D > 1.5, variations in relative water depth 


had no apparent effect on Cy or Cpe 


Recommendations for Further Research 


16. A phase difference may exist between the velocity 
of the ambient flow around a cylinder and the drag force. 
This suggests that values of Cy evaluated when the drag 
force was assumed to be zero may be in error. 

17. Further research was found necessary to determine 
and quantify the effects of vortex shedding on horizontal 
circular cylinders in oscillatory flows and to verify the 


relationship between C, and e/D for different Reynolds num- 


D 
bers and extended ranges of A/D. 
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APPENDIX A 


This appendix contains a computer routine for calcu- 
lating various wave parameters for experimental runs by the 
Airy wave theory. 


Nomenclature of results is as follows: 


HC = Wave height above cylinder (ft) 
HI = Wave height forward of cylinder (ft) 
T = Wave period (sec) 
LC = Wave length calculated by Airy wave theory 
(£t) 
L = Observed wave length (ft) 
H = Water depth (ft) 


E = Gap between cylinder and plane boundary (in) 
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VERTICAL FORCES AND PHASES 
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